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ABSTRACT 

We have developed a model for the polarization of curvature radiation for 
pulsars by taking into account the polar cap induced perturbation on the nonro- 
tating (slowly rotating) dipolar magnetic field, where the rotation effects such as 
aberration and retardation can be ignored. We have simulated a set of typical 
pulse profiles to understand the role of induced magnetic field on radio emis- 
sion of pulsars, and found to be significantly influencing the profile structure 
and polarization. Our model indicates that the intensity components and the 
polarization angle inflection point can get shifted to either leading or trailing 
side depending upon the prevailing conditions in the viewing geometry, the non- 
uniformity in source distribution (modulation) and the polar cap current induced 
perturbation. Also, we find an evidence for the origin of symmetric type circular 
polarization in addition to antisymmetric type. Our model predicts for a stronger 
trailing component compared to that on leading side of a given cone. 

Subject headings: polarization - pulsars: general - radiation mechanisms: 
non-thermal 



1. Introduction 

Pulsars are universally accepted as the fast rotating neutron stars with very strong 
magnetic field (~ 10^ — 10^^ G), which is predominantly dipolar. Their radio emission is 
thought to be curvature emission from relativistic plasma streaming out along the open field 
l ines of the dip o le those foot points on the neu t ron star surface define the polar cap region 



( ISturrockl Il97ll : iRuderman &: Sutherland! Il975l : ICheng fc RudermanI Il980l ). The simplistic 
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picture of pulsar radio emission that we can gather from the polar cap emission models is 
as follows. The rotationally induced very strong electric field component parallel to the 
magnetic field pulls out the primary plasma from the neutron star surface and accelerates 
them to very high energies of particles with Lorentz factor 7 ~ 10^ — 10^. These ultra- 
relativistic primary particles emit gamma-ray photons through curvature radiation which 
intern produce secondary electron-positron (e~, e^) pair plasma by the interaction with 
strong magnetic field. The secondary plasma further undergoes pair cascade process to 
develop into a high density plasma of 7 ~ 10^ — 10'^. Once the density of pair plasma become 
comparable to Goldreich- Ju lian charge density, it shields the parallel electric field component 
(iGoldreich fc Julian Ill969l ). Hence, the plasma motion become "force- free" along the field 
lines and it emits curvature radiation due to the curvature of the field lines. Since the 
brightness temperature of pulsar radio emission is very high (~ 10^^ — 10^*^ K), the emission 
is believed to be coherent emission from plasma bunches whose formation is believed to be 
due to plasma instability. 

Pulsars show high degree of polarization with systematic 'S'- shaped polarization po- 
sition angle (PPA) swing which is a characteri stic polarization property of p ulsar radio 
emission. The "rotating-vector-model" (RVM) of iRadhakrishnan fc Cookd (119691 ) attributes 
this characteristic 'S'- curve to underlying geometry, wherein the magnetic field is assumed 
to be mainly dipole field and relativistic beaming is in the direction of field line tangents. 
Ever since the highly successful RVM, fitting pulsar polarization profiles has been a standard 
model of constraining the underlying emission geometry such as the magnetic axis inclination 
angle with respect to rotation axis and the sight line impact angle with respect to magnetic 
axis. However, some pulsars show the polarization behavior that deviates from the standard 
'S' curve, partic ularly millisecond pu lsars, wherein the polarization sweep is noisy and on an 
average flatter (jXilouris et al.lll998f ). Several relativistic and plasma effects have been pro- 
posed to understan d these deviations which would modify t he PPA sweep through plasma 
propagation effects (IBarnard fc Aronslll986l : lMcKinnorull997l ). th rough distortion of underly - 
ing dipole field due to field aligned current flow above polar cap (IHibschman fc Aronsll2001[) 



through aberration of the beaming direct ion from strict parallelism (IBlaskiewicz et al 



1991 



Dyksl I2OO8I : iKumar fc Gangadhara II2012I ). t hrough distortion of under lying dipole field due 
to field aligned current flow above polar cap jHibschman fc Aronj2001), o r through multiple 
interacting orthogonal polarization modes ( McKinnon fc Stinebringj 1998 ). 



Phenomenologically, pulsar emission is recognized as central 'core' emission arising from 
the region near to the magnetic po le and surrou n ding ' cone' emission arising frorn con- 
centric rings around t he pole (e.g., iRankinI [ l983l. Il990l . Il993[ iMitra fc Deshpandel Il999l : 
Gangadhara &: Gupta! 1200 ll : iMitra &: RankinI |2002| ). However, due to an asymmetry in 



the phase location of components, there are some contrary arguments that the emission 
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is 'patchy' ( iLyne fc Manchester! Il988l ). Pulsars also show an asymmetry in the strength 
of the intensity components between the leading and trailing sides. The asymmetry ob- 



to the effects of rotation such as aberration and retardation, (e.g., Blaskiewicz e; 


al. 199ll: 


Hibsc 


nman & AronslboOl: Ganeadhara & Gup 


talboOl: Gu 


pta & Ganeadharall2003: 


Dvks et al. 


2004 


Thomas & Ganeadharall2007: Dvks.2008i: 


Dvks et al. 


2010: Thomas et allboiC 


: Kumar & Ganeadhar 


2012 


. Due to the rotationally induced aberration and asymmetry in the curvature of 



particle's trajectory, the inflection point of PPA profile (phase a t which PPA curve has 



maximum slope) lags the peak of the central intensity compoiient ([Blaskiewicz et al.l Il991 



Hibschman &: Arond 1200 ll : iDykd l2008l : iKumar &: Gangadhara 1 12012| ) , and th e leading side 

components become stronger and broader than the corresponding trailing o nes ( iThomas fc Gangadhara 
200?! : bvks et al.lboiol : iThomas et aPboiot iKumar fc Gangadhara I I2OI2I ). 



Han et al. 


1998; 


You & Han 


2006) 



importance from the point of view of emission mechanism. There are two types o f mod- 
els fo r the origin of circular po la rization: intrinsic to the emissio n mechanism (e.g.. 



1987 



Gil fc Snakowski 



1979 



1990al Jbl: iRadhakrishnan fc RankinI Il990l : iGangadhara 



Kuniar Sz Gangadhara 2OI2I ). and generated by the propagation effects (e.g., Cheng &: Ruderman 



1997 



Michel 
2010| : 



Melrosdl2003l ). Only the "antisymmetric"- type circular polarization (sign reversal un- 
der the component) was though t to be an intrinsic property of curvature radiation. Very 
recently iKumar fc Gangadhara I (|2012| ) have showed that the "symmetric"- type circular 
polarization (where the circular sign remains the same through out the component) is also 
possible within the framework of curvature radiation after taking into account of pulsar ro- 
tation and modulation (nonuniform plasma distribution). They have made an attempt to 
simulate the diverse behavior of polarization properties of pulsar radio emission by developing 
a relativistic model on full polarization of curvature radiation. 

Since the polarization behavior of pulsar emission is believed to depends upon the un- 
derlying geometry, any perturbation to the dipole field is expected to reflect in the po- 
larization profiles. The induced toroidal magnetic field due to field aligned, poloidal cur- 
rents above the polar cap pertur b the dipole field and it is quite significant in affecting the 
PPA (Hibschman &: Arond I2OOII ) and introduce phase shift into the intensity components 
Gangadhara ( 20051 ) . Due to polar cap currents, only an offset is introduced into PPA swing, 
but no phase shift in the inflection point. So far in the literature, only the PPA and the 
intensity components phase shift has been modeled. But the full polarization model that 
describes all the polarization parameters of radio emission in polar cap current perturbed 
dipole magnetic field do not exist. So, in an attempt to understand the polar cap current 
induced perturbation on the pulsar radio emission and polarization, we develop a complete 
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polarization model by taking into account of the polar cap currents. Our model considers 
the "force-free" motion of the relativistic plasma streaming out along the nonrotating per- 
turbed dipolar field lines emitting curvature radiation. We derive the viewing geometry for 
the perturbed dipole field in section 2. In section 3, we present the simulation of typical po- 
larization profiles for the cases of uniform and nonuniform distribution of radiation sources. 
In section 4 we give the discussion, and in section 5 the conclusion. 



2. GEOMETRY OF CURVATURE RADIATION IN THE PERTURBED 

DIPOLAR MAGNETOSPHERE 

In spherical polar coordinates (r, 9, (p) centered on the magnetic axis, the unperturbed 
dipole field is given by 

So ^ O) , (1, 

where n is the magnetic moment, 6 is the magnetic colatitude, and r is the radial distance 
from origin. Most of the polar cap models of pulsar emission include a current of relativistic 
charged particles streaming out along the open fiel d lines with intensity appro ximately equal 



to the Goldreich- Julian charge density (Jgj) (e.g.. iHibschman fc Aronsll200ll ). 

J = -^c;{n-Bo)Bo, (2) 
Ztt 

where (; = J/ Jqj is a scale factor of order unity, reflecting our ignorance of the actual current, 
and = VlVt is pulsar angular velocity. To simplify the things we treat <j to be a constant 
quantity through out this work. 

By assuming a small polar cap so that the variation in ■ Bq across the polar cap is 
also small, and that Bq is a point dipole with an axi symmetric current flow, th e induced 



magnetic field due to field aligned current is given by (IHibschman fc Arond 120011 ) 



„ / 2usin6' r \ , , 

-Bi= 0, 0, - ^ ^ cosa , 3 

V I' ^LC / 

which is purely toroidal and axisymmetric with respect to the magnetic axis, i.e, independent 
of the magnetic azimuth 0. The angle a is the magnetic axis inclination angle, and tlc = 
c/VL = cP/ (27r) is the light cylinder radius, where c is the speed of light, and P is the pulsar 
rotation period. Then the perturbed dipole field B = Bq + B^ with Bi <^ Bq, is also 
axisymmetric with respect to the magnetic axis. 



The differential equations for the magnetic field lines are 

dr r dO r sin 6 dO 

Br Bfj Ba, 



(4) 
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where Br, Bg, and are the respective r, 6, and (j) components of B in polar coordinates. 
Substituting the expressions for Br, Bg, and into Equation (jl]) yields the following three 
differential equations for the magnetic field lines: 

^ rfr = 2 cot ^ d9, (5) 



r 



sin 9 2q cos a 

and 



de = d(j), (6) 



1 'i"lc 

dr = d(j). (7) 



cos 9 q cos a 

Since Equation (JTj) can be deduced from Equations ([5]) and IQ, magnetic field hues can be 
defined uniquely by the later two equations. 

Integration of Equation (jS]) gives 

r = Ki sin^ 9, (8) 

where Ki is the integration constant. Using r = Vmax = re, the field line constant, at 9 = ir /2, 
we get Ki = Te- Substituting for Ki in Equation (jH]) gives 

r = Te sin^ 9. (9) 

Substitution of Equation ([9]) into Equation (E]) and integration gives 

re cos ^ = + 7^2, (10) 

l<; cos a 

where K2 is the integration constant. Using (j) = (pi aX 9 = 9i, we get K2 = cos 9i — 
tlc 0i/(2<;cosa), where 9i and 0^ are initial colatitude and azimuth, respectively. Substitut- 
ing for K2 in Equation ( ITOll yields 

2re 

= 0j H ^cosa (cos6' — cos6'j). (11) 

The position vector r of any point alon g the perturb e d dip olar field line is given by the 



general expression given in Equation (1) of iGangadharal (120101 ). where the parameter has 



to be replaced by the expression that we have given in Equation (fTTI) . By using a = 30°, 
magnetic axis rotation phase 0^ = 30°, re = 1 rLc, -P = 1 s, <^ = 1, 6*, = 0°, 0j from 0° to 360° 
in steps of 45°, and r from to 0.6 tlc we have traced the field lines and are given in Figure 
[H Dashed line curves represent the field lines of an actual dipole while the solid line curves 
represent those of the perturbed dipole. Polar cap currents induce the toroidal component 
of magnetic field and hence there is an azimuthal twist of the field lines of the perturbed 
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field with respect to the magnetic axis. Note that the field fines of tfie unperturbed field 
have no azimuthal twist as the field is purely polar with respect to the magnetic axis. 

The curvature of the perturbed dipolar field line is given by 

db 1 db , . 

^ Ts ^ Wide' 

where ds — \h\d9 is the arc length of the field fine, b — b/\b\, and b — dr/dO is the tangent 
vector to the field fine. The magnitude of b can be readily found as 

|6| = ^^^^ VlO + 6 cos 29 + 4 A2 sin^ 0, (13) 

where the parameter A = 2^r cosa/rLc- 
Then, the radius of curvature p — l/|fc| of field line becomes 

p= ^, ^-""''^ (14) 

where k = db/ds is the curvature vector of the field line, 

Di = 2 (9 + ISA' + A^) + (5 + A') d^, 
L>2 = 24:A^dl+ {3- A^)dl 
Ds = 6(1-A2)M2, 



and 

di = 10 + 6cos2^ + 4A2sin2e, 
d2 = 6(1- A^) sin 2^. 

The velocity of rclativistic source (particle or plasma bunch) which is constrained to move 
along the perturbed field line is given by 

V — K cb, (15) 

where the parameter k = \J\ — 'Vj'i^ with 7 being the source Lorentz factor. 

Then the acceleration a — dv/dtol 'force- freely' streaming plasma due to curvature of 
field lines is given by 

a = {Kcfk, (16) 

where we have used ds — \h\d6 — kc dt. 
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The accelerating relativistic particles emit beamed emissions in the direction of their 
velocities. At any given rotation phase (j)' = 0^ and emission altitude r, observer receives the 
emissions from a finite beaming region constituting the flux of magnetic field lines. Observer 
receives maximum beamed emission when sight line fi exactly aligns with source velocity v. 
Observer sight line n = (sin^, 0, cosC), where the angle ( = a + a with a being the sight 
line impact angle with respect to magnetic axis. Since the beaming angle of the beamed 
emission from the relativistic particle is ~ I/7 about its velocity, beaming region boundary 
can be specified by the emission points at which n ■ b = cos(l/7). Therefore the emission 
points within the beaming region satisfy the condition — l/7<?7<l/7, where t] is the angle 
between n and b. Following Kumar & Gangadhara (2012), we define the emission point 
coordinates: magnetic colatitude and azimuth of the emission points within the beaming 
region as 6^ and (pe, respectively, wherein at the center of the beaming region Of. = 6q and 

4>e = 00- 

Next, for a given emission altitude r, we find the emission point coordinates ^0 and 
00 of the beaming region center and the range of 6 and which define the beaming region 
boundary. Solving rh ■ b = cos r = (1 + 3 cos 26) / VlO + 6 cos 29 + 4A2 sin^ 6 for 6 gives 

cos^r + cosr^J^l - cos^ r + 8 (^1 + - 1 j , (17) 

where r is the angle between m and b. But the angle F between rh and n is given by 

cos F = cos a cos C + sin a sin ( cos 0'. 

For an exact alignment of h and b at the beaming region center, h ■ b = 1, which implies 
r = F. Hence we have 




cos 29o = ^y^l- cos^ F + cos Fy (^1 - cos2 F + 8 (^1 + ^ 

By solving n x 6 = 0, we obtain 

, cscF(sinC(Acosacos0' — 3cos6'sin0') — AsinacosO , . 

sm0o = (19) 

VA2 + 9 cos2 e ^ ' 

and 

cscF(sinC(3cosa;cos6'cos0' + Asin0') — SsinacosCcos^^) 

cos 00 = . (20) 

VA2 + 9cos2 ^ ' 

In Figure |2] we have plotted the emission point coordinates 6q and 0o of the beaming re- 
gion center and corresponding radius of curvature p of source trajectory at those points as 
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functions of rotation phase cj)'^ of the magnetic axis by choosing the parameters a = 30°, 
cr = ±5°, emission altitude r„ = r/ri^c = 0.1, P = 1 s, and = 1. The dotted hne curves are 
for the parameters in the case of unperturbed dipole and the sohd hne curves represent those 
in the perturbed case. The colatitude is found to be not noticeably get affected whereas 
there is an upward vertical shift in azimuth 0o due to the induced toroidal magnetic field. 
Therefore at any 0^ radiation from the emission points aberrated in the direction opposite 
to the direction of pulsar rotation for positive a (i.e., for ( > a), where as for negative a 
(i.e., for ( < a) it is aberrated in the direction of pulsar rotation. In other words emission 
locations in get shifted to later phases with respect to those due to unperturbed dipole for 
the positive a, whereas for negative a, they get shifted to earlier phases (the phase shift of 
an antisymmetric point of 00 with respect to (p'^ = is indicated by longer arrows). Note 
that the inflection point of 0o (phase at which 0o curve has maximum slope) is found to 
be remains unaffected as indicated by shorter arrows. Since the effect of induced toroidal 
component of the magnetic field due to polar cap currents is just a twist of the field hues 
around magnetic axis and same on leading and trailing side field lines, p in the perturbed 
configuration will remain almost symmetric with respect to 0^ = 0° similar to that in the 
unperturbed one but with smaller values because of an induced curvature. 

Next at = 00, we solve n -b = cosr]max = cosl/7 and find the allowed range of r: 
F — l/7<r<r + l/7, and hence deduce the allowed range of 6 from Equation f|T71) . Then 
for any 6 within its allowed range, we find the range of 0: 0o — 50 < < 0o + S(t> by solving 



n-b = cos 1/7. Note that the expression obtained for 6(f) is same as that given in lGangadhara 



(I2OIOI ) except for the fact that the angle r has to be replaced with the one that we have given 
above. We noticed that the behaviors of beaming regions with respect to rotation phase 0^ 
are almost similar to those in the case of unperturbed dipole of non-rotating magnetosphere 



(lGangadharall2010l ) . Note that the whole beaming regions of the perturbed dipole get shifted 



in the coordinate with respect to those of the unperturbed dipole (see Figure [2]). 



3. POLARIZATION STATE OF THE RADIATION FIELD 



The spectral distribution of the radiation field emitted by the accelerating relativistic 
plasma at the observation point Q is given by (Jackson 19981 ) 



1 qe 



iLuRo/c 



271 Ro C 



\h\ hx[{n- (3) X /3] 



KC 



(21) 



Note that time t in the above equation has to be replaced by the expression given in Equation 
(6) of Kumar & Gangadhara (2012). Following the method given in Kumar & Gangadhara 
(2012) we solve the integral of Equation ( l2T]) and estimate the polarization state of the 
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radiation fie l d in t erms of tlie Stokes parameters /, Q, U, and V, wliicli are defined in 



Gangadharal (120101 ) 



3.1. Emission With Uniform Radiation Source Density 

By assuming uniform radiation source distribution in tlie emission region of magneto- 
spfiere we estimated the polarization state of the emitted radiation. The simulated contour 
plots of total intensity /, linear polarization L = y^Q^~+TP , and circular polarization V 
in the perturbed dipole configuration are given in Figures |31 The parameters used for sim- 
ulation are a = 30°, a = ±5°, r„ = 0.1, P=ls, <^ = 1, 7 = 400, u = 600 MHz and 
rotation phase 0^ = 0°. Note that the parameters in each panel are normalized with the 
corresponding peak value of total intensity. Also note that the contour plots are plotted in 
such a way that the emissions from the leading side to trailing side of the beaming regions 
in azimuth fall over the region from left to right side of panels. 

The polarized emissions from the beaming region gets rotated as indicated by the con- 
tours patterns in [6, (j))— plane where as those due to unperturbed dipole will be either 
symmetric (/ and L) or antisymmetric {V) with respect to {6, (po)— plane. The polarized 
state of the emission from the beaming region in the unperturbed dipole geometry is given 
in Kumar & Gangadhara (2010). The rotation of contour patterns is due the perturbation 
to the actual dipole geometry due to polar cap currents. Also it introduces an asymmetry 
into the beaming region emissions in such a way that there is a selective enhancement of 
emissions over larger 6 and larger region due to a larger resultant curvature. This is very 
clear from the V parameter contour plots where there is a selective enhancement of positive 
circular over the negative circular. Note that the direction of rotation of contour patterns 
in {9, (j))— plane is opposite between the positive and negative a cases. It is in opposite 
direction to the one with rotating unperturbed dipole for positive a while for negative a 
it is in same direction (see Kumar & Gangadhara 2012). Also note that the magnitude of 
rotation of intensity patterns is found to be almost independent of rotation phase of the 
magnetic axis unlike the one due to effect of pulsar rotation where it noticeably decreases at 
outer phases. 

The net emission that the observer receives at any instant of time will be an incoherent 
sum of emissions from diff erent plasma particles/bunches within the beaming region. Fol- 



lowing [G^ji^adharal ( 120101 ). we have simulated the polarization profiles: total intensity Ig 



linear polarization Lg = ^/Ql + U^, and circular polarization Vs, and polarization position 
angle ips = (1/2) tan~^(?7s/Qs) and plotted in Figure HI For the simulation, we have used 
the parameters the same as in Figure [3] and assumed that the source density is uniform 
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through out the pulsar emission region. The emission on leading side (0^ < 0°) is roughly 
same as that on trailing side (0^ > 0°) due to the symmetry of curvature of field lines on 
both the sides (see Figure [2]) which is same as in the unperturbed dipole geometry. Even 
though emissions from individual sources are highly polarized, the net emission from the 
whole beaming region becomes less polarized due to incoherent addition. However a small 
quantity of positive circular polarization survives due to the asymmetry in the positive and 
negative polarities (see Figure [3]) which would become more important when the source den- 
sity becomes nonuniform in the emission region. The polarization position angle (ips) curves 
(solid line curves) shift upward with respect to those due to unperturbed RVM model (dot- 
ted ones) o nly, and the PPA inflection point almost left unchanged. For comparison, PPA 
derived by iHibschman &: Arond ( 120011 ) are also superposed wit h our simulations arid both 
are in good agreement (see the dashed line curves). Note that IHibschman fc AronsI (l200l[ ) 
considered the emission from the beaming region center only, and the simulations show that 
the contributions from neighboring emission points will have a least effect on net PPA. How- 
ever they become important once there is a density gradient in the plasma distribution as 
we show later. 



3.2. Radiation Emission from Nonuniform Distribution of Sources 

The sub-pulses in pulsar individual pulses suggest that the emission region of pulsar 
magnetosphere may not be filled with uniform plasma density. When the observer sight line 
encounters plasma columns that associated with the sparks activity on polar cap, components 
like features in intensity can result. Hence the net emission that the observer receives will 
be modulated emission over the uniform background emission. We model the nonuniform 
source distribution in the emission region in terms of "modulation" function as considered 
in Kumar & Gangadhara (2012). The intensity components of average profiles are assumed 
to be nearly Gausians in shape, and hence they have been fitted with appropriate Gaussians 
(e.g., Kramer et al.lll994[ ). We define the modulation function as 

fiOA) = fofeU, (22) 

where /o is the amplitude, and fg = exp [—{9 — Opj^/ag] and = exp [—{(f) — 0p)^/c'"^] are 
the Gaussian modulations in 6 and (p, respectively. The parameters [dp, (pp) define the peak 
location of the Gaussian and ae = wg/{2\/]n2) and = w^/ (2A/ln2), where wg and w^ 
are the corresponding full width at half-maximum (FWHM) of the Gaussian in the two 
directions. 
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3.2.1. Emission With Azimuthal Modulation 

By considering a modulation in the magnetic azimuthal direction having peak location 
in the magnetic meridional plane, in the two cases: (f)p = 0° for a = 5°, and (pp = 180° 
for a = —5°, we have computed the polarization profiles of observed radiation and given in 
Figures & El respectively. For simulation we have used the parameters /o = I, fg = 1, and 
other parameters the same as of Figure |3l In the case of cr = 5° (see Figure |5]) intensity peak 
shifted to later phase with respect to 0^ = 0°, whereas for cr = — 5° it gets shifted to earlier 
phase (see Figure [6]). 

This is because, due to the effect of induced magnetic field by the polar cap currents, the 
dipolar field lines get perturbed. As a result, the emission locations are shifted in magnetic 
azimuth as per Figure |2l and hence intensity get phase shifted to the later and earlier phases 
for the positive and negative sight lines, respectively. The phase shifts of the intensity 
peaks in the case of a = 5°, and for cr^ = 0.1 and 0.4 are found to be 0°.54 and 7°. 22, 
respectively, whereas those in the case of a = —5°, and for = 0.1 and 0.25 are found to be 
0°.75 and 1°.50, respectively. The increase in phase shift of intensity peaks with increasing 
can be explained as follows. In the perturbed dipole geometry the emissions locations 
(particularly in magnetic azimuth) get phase shifted with respect to the unperturbed one. 
Because of this phase shift, the curvature of the source trajectory becomes asymmetric about 
the antisymmetric points of the magnetic azimuth (i.e., 0o = 0° or 180° phase locations for 
positive and negative a respectively). For cr = 5° curvature becomes larger on the trailing side 
of 00 = 0° compared to that on the leading side, whereas it is vice versa about 0o = 180° 
for a = —5° (see Figure [2]). Therefore the phase shifts of the modulated intensity peaks 
will be larger as compared to those of the antisymmetric points of 0o (where the modulation 
peaks are located) and this extra phase shift becomes larger for broader modulations. Further 
intensity component gets broader as well as shifted more for the negative a as the modulation 
mapped on to broader region in observer frame compared to that for the positive a (see the 
case of 0",^ = 0.1 of Figures [5] & [6]) . 

Circular polarization Vg is antisymmetric for the steeper modulation (cr^ = 0.1), and 
changes sign from negative to positive in the case of cr = 5° and from positive to negative 
in the case of cr = —5°. In both the cases of cr = 5°, and —5° positive circular becomes 
stronger compared to the negative circular due to an enhancement of the unmodulated 
positive circular over the negative circular within the beaming region (see Figure [3D. Also 
the phase location of the sign reversal of circular either leads or lags the intensity peak due 
to above said asymmetry. In the case of broader modulation (cases = 0.4, and 0.25 of 
Figures [5] and [6l respectively), circular becomes almost positive through out the pulse. Hence 
it becomes almost symmetric type as the effect of modulation become not so significant over 
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the beaming regions in comparison to the induced asymmetry between the unmodulated 
positive and negative circulars that resulted from polar cap currents. 

Linear polarization almost follow the total intensity Is profile with lower magnitude 
due to incoherent addition of emissions within the beaming region. Further, it is slightly 
enhanced on leading side, where circular is smaller, than on trailing side in the case of cr = 5°, 
and vice versa in the case of o" = —5°. PPA is increasing PPA( "counter clockwise" or 
"ccw" swing as dips/dcj)'^ > 0) in the case of a = 5°, and decreasing PPA ('clockwise" or "cw" 
swing as dips/d(j)'^ < 0). Hence the correlation of the sense reversal of circular polarization 
from negative to positive with the increasing PPA and vice versa, in the case of pulsars with 
antisymmetric type circular polarization. Inflection point of PPA (phase at which dips/ d(f)'^ 
is maximum) is not noticeably get affected for the broader modulations (see cases = 0.4 
and 0.25 of Figures E] and O respectively), similar to the one with an unmodulated emission 
(see Figure Hj). But for steeper modulation (cases = 0.1 of Figures [5] and [6]) , it is shifted 
to earlier phase in the case of cr = 5° and to later phase in the case of a = —5° which is 
opposite to the intensity phase shifts. Note that the phase shift of PPA inflection is found 
to be shifted towards the stronger side. The phase shift of PPA inflection point in the 
cases of (T = 5° and —5° with a^j, = 0.1 are found to be 0°.73 and 0°.63, respectively. The 
smaller phase shift of PPA inflection point in the case of a = —5° as compared to that in 
the case of a = 5° is due to effectively broader modulation in the observer frame, and hence 
smaller effect of modulation over the beaming regions emission in the case of cr = — 5°. 

For the next case we consider the modulation symmetrically located on either sides with 
respect to meridional plane: (f)p = (180 ± 30)° for a = —5° and 0p = ±40° for a = 5°. The 
simulated polarization profiles are given in Figure [3 The modulation strength / that the 
observer finds at the beaming region center and the corresponding radius of curvature p of 
particle's trajectory are also given. Even though p is roughly symmetric about 0^ = 0, 
it becomes asymmetric in the strength at the modulation peaks between the leading and 
trailing sides. The p at the leading side modulation peak becomes smaller than that on the 
trailing side in the case of a = —5°, and vice versa in the case of cr = 5°. This is due to the 
phase shift of the modulation peaks (mid point of the modulation peaks is indicated by the 
thick arrows) due to the phase shift of the emission points (see Figure [2]). As a result an 
asymmetry arises in the strength of intensity components between leading and trailing sides: 
component become stronger on leading side than the trailing one in the case of a = —5° and 
vice versa in the case of cr = 5°. Also, since the p at the leading side modulation peak becomes 
less steeper than that on the trailing side in the case of cr = —5°, the leading component 
becomes broader, and vice versa in the case of a = 5° as indicated by the second row panels. 
In both the cases of cr, similar to the modulation the outer phase emissions become stronger 
than the inner ones, and K clearly shows that. However the larger asymmetry between the 
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positive and negative circulars under the leading side component in a = —5° case and that 
under the trailing one in o" = 5° case is due to the selective enhancement of positive circular 
of the unmodulated emissions (see Figure [3]). 

3.2.2. Emission With Polar Modulation 

By assuming a hallow cone modulation around the magnetic axis, i.e, plasma density 
gradient only in the polar direction, we have simulated the polarization profiles and given 
in Figures E] and O The parameters used for the simulation are /q = 1, = 1, and the 
rest are the same as of Figure [31 In the cases of Op = 3° and ag = 0.005, and 0.002, since 
the minimum of ^ ~ {2/3)a = 3°3 the observer sight line just grazes the hallow emission 
cone and hence a single intensity component is observed. But in the case of = 4° and 
ag = 0.003 observer cuts the hallow cone twice and hence results in double components. In 
all cases, intensity peak (/^ peak in the cases of 9p = 3° and ag = 0.005, and 0.002 and mid 
point of the intensity peaks in the case of 6p = 4° and ag = 0.003) is more or less remains 
at 0^ = 0°. This is because, even with the perturbation of the dipole, both the colatitude 
and the corresponding radius of curvature remain symmetric about (p'^ = 0°. However there 
is a small phase shift of PPA inflection point in the cases of 6p = 3° and ag = 0.002, and 
6p = 4° and ag = 0.003, but the shift is in the opposite direction between a = 5°, and —5° 
cases. 

In the case of 6p = 3° and ag = 0.005, a tiny circular survives due to not so steeper 
modulation, but in the cases of 6p = 3° and ag = 0.002, and 6p = 4° and ag = 0.003 it 
is quite significant. In the case of 6p = 3° and ag = 0.002 where the observer sight line 
grazes the modulation at larger colatitude, a symmetric type negative circular observed due 
to the selective enhancement of beaming region emissions over the smaller 6 part (see Figure 
|3]). But in the case of 9p = 4° and ag = 0.003, an antisymmetric type circular observed 
with sign changing from negative to positive under leading component and vice versa on 
trailing. This is due the selective enhancement of emissions over either smaller values of 6 
or larger depending upon when the sight line crosses the hallow cone. Note that with the 
modulation in 6, the sign reversal of circular is same for both the cases of a unlike the case 
with modulation in (see Figures |5] and |6]) because of the similar shape of the colatitude 
in both the cases (see Figure [2]) . Since the beaming region emission due to uniform source 
distribution have larger positive circular than the negative circular (see Figure [3]), modulated 
profiles also show slightly larger positive circular than the negative circular. 
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3.2.3. Emission With Modulation in both the Magnetic Colatitude and Azimuth 

In general, the plasma distribution can be nonuniform in both the polar and azimuthal 
directions and hence we consider the modulation in both the 6 and 0. The extreme cases of 
modulation: modulation in (p dominating over that in 9 for <^ ag and the modulation in 6 
dominating over that in (p for ag <^ cr^ are detailed in sections 3.2.1. and 3.2.2. respectively. 
The simulations with the combinations of modulation parameters used in the sections 3.2.1. 
and 3.2.2. give the general polarization profiles. 

For the case with Gaussian modulation with peaks located symmetrically on either sides 
of the magnetic meridional plane in the conal ring of colatitude 4°, we consider two sets of 
azimuthal peak locations at </)p = ±50°, and ±30°. The simulated polarization profiles are 
given in Figure [TO] for the case of o" = 5°. For simulation we used /o = 1, cr^ = 0.01, = 0.1, 
and the other parameters the same as of Figure [31 Since the phase shift of emission points 
due to the effect of polar cap currents is only in (p, but not in 6, the observer sight line 
selectively encounters the modulation depending upon their azimuthal locations as shown 
in Figure [TOl The sight line encounters stronger modulation on leading side as it passes 
closer to the Gaussian peak than that on the trailing side in the case 4>p = ±50°, and vice 
versa in the case (pp = ±30°. However, the asymmetry in the strength of the components 
between the leading and trailing sides become enhanced or reduced in the intensity profiles 
depending upon the asymmetry in the curvature at the modulation peaks as explained earlier 
for Figure [71 The p at the modulation peak on the leading side become larger than that on 
the trailing side. Hence the decrease in the asymmetry in the strengths of the modulated 
intensity components between the leading and trailing sides in the case (pp = ±50°, where as 
it increases in the case (j)p = ±30°. The inner circulars become stronger than the outer ones 
in the case (pp = ±50° as the observer selectively encounters the modulated regions which 
are closer to the magnetic meridional plane, and vice versa in the case 0p = ±30°. Further, 
it is more pronounced on the trailing side than the leading side due to an asymmetry in the 
modulation encountered by the observer. 

Note that if one considers modulation in 6 direction which is even steeper, then the 
circular becomes almost symmetric type with only the inner circulars sign for (f)p = ±50°, 
and only the outer circulars sign for (j)p = ±30°. Also note that if one chooses negative 
sight line a = —5° and the modulation peaks situated symmetrically in azimuthal direction 
with respect to the meridional plane, then the behavior of the strength of modulation that 
the sight line encounters and intensity components on the leading and trailing sides will be 
opposite to the case of a = 5°. 
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4. Discussion 

The distortion of magnetic dipole field due to the field aligned polar cap currents have 
a significant role on the pulsar emission and polarization. The induced toroidal magnetic 
field due to the poloidal currents above the pulsar polar cap perturbs the underlying dipole 
field. The resultant field is also an axisymmetric with respect to magnetic axis wherein 
the field lines are twisted around the magnetic axis due to an azimuthal component. Since 
the plasma particles (or bunches) are constrained to move along the field lines with which 
they are associated, the beaming direction of the relativistic emission gets aberrated in the 
direction of induced twist of the field lines. As a result, for the case of positive sight line 
(o" = +ve), the emission points of the perturbed dipole field are shifted to later phases (or 
later times) in the magnetic azimuth than those of the unperturbed dipole field, and vice 
versa for the case of negative sight line {a = — ve). Since the induced twist of the field lines 
is around the magnetic axis only, the emission points are found to be more or less unaffected 
in the magnetic colatitude. As the perturbed dipole field remains axisymmetric, for a given 
field line constant r^, the radius of curvature of the perturbed dipole field lines become 
function of colatitude only. Hence the curvature remains symmetric about 0^ = 0, which 
is similar to the magnetic colatitude. Therefore the emission due to the uniform radiation 
source distribution and its polarization profiles remain symmetric about 0^ = as shown 
in Figure HJ which is similar to the emission in unperturbed dipole field except the survival 
of a small circular polarization due to polar cap current induced asymmetry. 

However, in emissions due to the nonuniform source distribution with a density gradient 
effectively in the magnetic azimuth, the modulated intensity components get shifted in phase. 
Since the phase shift of the magnetic azimuth between the two cases [a = ±) are of opposite 
directions, the modulated intensity components shifts to the later phase in the case of positive 
a and to the earlier phase in the case of negative a. Even though the radius of curvature 
p is symmetric about 0^ = 0, it becomes asymmetric about the shifted modulation peak 
where the observer encounters the maximum plasma density. Hence the phase shift of the 
modulated intensity components become significantly different from that of the peak location 
of the broader modulation. But the intensity components do not get shifted in phase if the 
modulation is effective only in the polar direction as both the 6 and p remain symmetric 
about 0^ = 0. 



Hibschman &: AronsI (1200 ll ) showed that due to the distortion of the underlying dipole 
field, the resultant PPA curve gets upward shift in the PPA v/s rotation phase diagram 
compared to the one with the unperturbed dipole field, leav ing its inflection poin t rema ins 



roughly unchanged. Our simulations confirm the results of iHibschman fc AronsI (120011 ) in 



the case of unmodulated emission. However, once the radiation source distribution becomes 
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nonuniform with larger density gradient, as we showed in sections 3.2., the phase shift of the 
PPA inflection point becomes quite significant. Note that, if there are phase shifts between 
the intensity components and that of the PPA inflection point then they will in the opposite 
directions. 

Due to the phase shift of the emission points in the perturbed dipole field in the mag- 
netic azimuth but not in the colatitude, the observer, in general, encounters an asymmetry in 
the strength of the two directional (in both the polar and azimuthal directions) modulation 
between the leading and trailing sides as shown in Figure [HH Depending upon the azimuthal 
peak location of the Gaussian modulation in a given cone, the sight line selectively encoun- 
ters the modulation either on leading side or trailing side. Hence, an asymmetry arises in 
the strength of intensity components between the leading and trailing sides. However, the 
asymmetry in the intensity profiles could get enhanced or reduced compared to that of the 
modulation. This depends upon the asymmetry in the curvature of source trajectory at 
the modulation peaks between the two sides which is induced due to the phase shift of the 
modulation peaks. Note that, even with modulation only in the azimuthal direction, an 
asymmetry arises in the strength of intensity components between the leading and trailing 
sides as shown in Figure [3 

Statistically, even though the leading side components dominate over the trailing ones, 

the re verse is also quite significant. Of the 88 cone dominated pulsars studied by lLyne fc Manchester 
( 119881 ) . 53 of the pulsars clearly have a stronger leading component th an the trailing one, 
and the 24 remaining show the reverse behavior (IBlaskiewicz et al.lll99l[ ). Our pulse profiles 
simulation also shows that depending upon the emission beam geometry and the modulation, 
the leading side component could become either stronger or weaker than the corresponding 
trailing one. However, in the literature, the asymmetry between the leading and trailing side 
components strength has been attributed to the effects of rotation. The weakening of the 
trailing side component is believed to be caus ed by the smaller curvature of source trajectories 



due to rotationally induced asymmetry (e.g . . IBlaskiewicz et al 



1991 



Thomas fc Gangadhara 



20071 : iDvks et al.ll2010l : iThomas et al.ll2010l : Kumar fc Gangadhara II2012I ): and the erihance - 
ment of trailing side emission due to squeezing into a narrower component ( iDyks et al.ll2010l ). 
However, the polar cap currents combined with pulsar rotation could be responsible for the 
asymmetries found in components strength. 

Also, due to the polar cap current induced asymmetries in the emission points phase shift 
and the curvature of the source trajectory in the observer's frame, an asymmetry introduced 
in the width of the intensity components, as shown in Figure [71 The components on leading 
side could become either broader or narrower than those on the trailing side depending upon 
the viewing geometry, the polar cap current induced asymmetries, and modulation. Statisti- 
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cally, of the 24 pulsars studied by lAhmadi fc Gangadharal ( 120021 ) . 19 pulsars have a broader 
leading com ponent compared to tha t on t h e trailing si d e. By taki ng into account of pul- 



sar rotation 



Ahmadi &: Gangadharal (120021 ) , iDyks et al.l (120101 ) and iKumar fc Gangadhara 



( 120121 ) have confirmed that the leading side components becomes broader than the corre- 
sponding trailing ones of the same cone emission due to rotationally induced asymmetry in 
the curvature of source trajectory between the two sides. However, the reverse case might 
be due to polar cap current induced effects. 

We showed that both the "antisymmetric" (with an asymmetry in the strength between 
the opposite polarities, in general) and "symmetric" types of circular polarization could be 
observed any-where within the pulse window. This is possible due to an induced asymmetry 
in the polarized emission from the beaming region, such as asymmetry between the positive 
and negative polarities, and asymmetry due to the rotation of circular contour patterns in 
(6*, 0)— plane, and the modulation of emission due to nonuniform distribution of sources. 

Our simulations show that the sign reversal of circular polarization near the central 
region of the pulse profile from negative to positive is correlated with the increasing PPA 
(wherein the PPA curve has positive slop e) and vice versa. This beha vior is found in mari y 
puls ars (IRadhakrishnan fc Rankin! Il990l ) , and has been confirmed by iGangadharal ( 120101 ) , 



and lKumar fc Gangadhara 



( 120121 ). However, there is also an association of the sign reversal 



of circular polar i zation from posit ive to negative with the increasing PPA, and vice versa 



( iHan et al 



19981 : IYou fc Hani |2006| ) . In the case of pul sars with the sign r eversal of circular 



polarization not associated with central core region, IYou fc Hani ( l2006l ) have found that 
both the sign reversals of circular from negative to positive and from positive to negative 
are associated with both the increasing or decreasing PPA, and our simulation of double 
component pulse profiles confirms this behavior. In pulsars with the "symmetric" type of 
circular polarization, we observed no c orrelation betwee n the sense of circular polarization 
and the PPA swing. Hence we confirm iHan et al.l ( 1l998l ) findings that the positive circular 
could be associa ted with either increasing or de creasing PPA, and similarly negative circular 
too. Note that iKumar fc Gangadhara I (120121 ) have also confirmed these findings for the 
emission in the rotating dipole field. 

In this work we have considered the altitude independent polar cap current density, 
but in reality it would have an altitude dependency due to the diverging nature of dipolar 
open field lines and conservation of escaping plasma. Since we have considered constant 
emission altitude r across the pulse profile, our simulations are still valid in explaining the 
trend of the effect of polar cap current on pulsar polarization profiles. However, if one 
consi ders the emissions frora the range of altitude that increases from inner to ou ter phases 
(e.g., iGangadhara &: Guptal I2OOII : iGupta &: Gangadharal 120031 iDyks et al.l 12004 ). then the 
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altitude dependent current density is more appropriate one. Since the understanding of polar 
cap current induced effects on pulsar emission and polarization in the nonrotating perturbed 
dipole field is needed for clarity, we have not considered the pulsar rotation here. However, 
the combined role of polar cap currents and pulsar rotation on the pulsar radio emission and 
polarization has to be investigated. 

5. Conclusion 

By considering the perturbation of the dipole field due to the field aligned polar cap 
currents we have developed a model for pulsar emission and polarization. By incorporating 
the nonuniform distribution of radiation sources in the emission region of the distorted dipole 
field, we have simulated a few typical polarization pulse profiles. Based on our simulations, 
we conclude the following points: 

1. Due to the polar cap current induced perturbation of the underlying dipole field, the 
emission points in magnetic azimuth get shifted to later phases in the case of positive 
(T, and to earlier phases in the case of negative o compared to those due to unperturbed 
one (see Figure [2]). 

2. Intensity components those resulted due the modulation that dominates in azimuthal 
direction over that in polar direction shift to either later phases or earlier depending 
upon the viewing geometry. The magnitude of phase shift depends upon the steepness 
of modulation and the polar cap current induced asymmetry in the curvature of source 
trajectory about the modulation maximum. 

3. Due to the polar cap current induced asymmetry and modulation, there arises a small 
phase shift in the PPA inflection point. 

4. The leading side intensity components could become stronger and broader than the 
corresponding trailing ones or the vice versa depending upon the viewing geometry 
and azimuthal location of the Gaussian modulation peaks in the concentric cones. 

5. Both the "antisymmetric" type circular polarization with an usual asymmetry in the 
strength of the opposite polarities, and the "symmetric" type circular polarization 
become possible anywhere within the pulse window. It depends upon the viewing 
geometry and modulation in the emission region of the perturbed dipole field. 
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Fig. 1. — Geometry of curvature radiation from a relativistic source S in a distorted (or 
perturbed) magnetic dipole due to polar cap currents. Magnetic axis rh is inclined by an 
angle a with respect to pulsar rotation axis Q. The dashed line curves represent the field 
lines of a pure dipole whereas the solid line curves represent those of the perturbed dipole. 
The observer's sight line h lies in the fiducial plane (XZ-plane), and the observation point Q 



is at a distance R from the emission point Pg. We have used a — 30°, 0' 



30° 



1 r^c, 



P = 1 s, ^ = 1, 0j from 0° to 360° in steps of 45°, Oi = 0°, and r from to 0.6 tlc to sketch 
the field fines of the pure dipole and those of the perturbed dipole. 
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Fig. 2. — Magnetic colatitude and azimuth 0o of the beaming region center are plotted 
as functions of rotation phase (p'^. The radius of curvature p is plotted in lower panels. The 
dotted and solid line curves are for the unperturbed and perturbed cases of dipolar magnetic 
field. The parameters used are r„ = r/ri^c = 0.1, P = 1 s, and ^ = 1. In the 0o panels, the 
long and short arrows indicate the antisymmetric and inflection points of (po in the perturbed 
dipole, respectively. 
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Fig. 3. — Polarization state of the emissions from the beaming regions with uniform radiation 
source density at the rotation phase 0^ = 0°. For simulation we have used the parameters 
a = 30°, a = ±5°, r„ = 0.1, P = ls, ^ = 1, 7 = 400, and v = 600 MHz. 
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Fig. 4. — Simulated polarization profiles for emission with uniform source distribution. In 
the PPA (tps) panels, dotted line curves represent the RVM for the unperturbed dipole, 
wher eas dashed and solid o nes re present the PPA in the perturbed dipole geometry derived 
from iHibschman &: AronsI ( 120011 ) and our calculations, respectively. The solid arrows point 
to the PPA inflection points from our simulation. The parameters chosen for simulation are 
same as in Figure |3l 
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Fig. 5. — Simulated polarization profiles for the emission from the nonuniform distribution 
of sources (modulation) in azimuthal direction for a = 5° case. The parameters used for 
simulation are /o = 1, /e = 1, (pp = 0°, and other parameters are the same as in Figure O 
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a=30°, o-=-5° 




Qr=30°, o-=5° 





-20 -10 10 20 _i5_i0 -5 5 10 15 



Fig. 7.— Same as Figure E] except 0p = (180 ± 30)° for a = -5° and 0p = ±40° for 
o" = 5°. In the top row panels, the thick hne curve represents the modulation strength / as 
seen by the observer at the central point of the beaming regions whereas the thin line curve 
represents the corresponding radius of curvature p; the thick line arrows mark the midpoint 
of modulation peaks whereas the thin line arrows mark the symmetric point of p. Both / 
and p are normalized with their corresponding maximum values. 
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Fig. 8. — Simulated polarization profiles for emission with modulation in polar direction for 
the case of cr = 5°. Here /o = 1, /</, = 1, and the other parameters are the same as of Figure 

El 




^'"("^ r„n r^n 



Fig. 9. — Same as Figure [8] except for cr = —5° and 0p = 180°. 



- 29 - 




Fig. 10. — Simulated polarization profiles for emission with modulation in both the polar 
and azimuthal directions for the case of ex = 5°. Here /o = I, 9p = 4°, = 0.01, = 0.1, 
and the other parameters are the same as of Figure [31 



